Optoelectronic Properties and Excitons in Hybridized 
Boron Nitride and Graphene Hexagonal Monolayers 
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We explain the nature of the electronic band gap and optical absorption spectrum of Carbon 
- Boron Nitride (CBN) hybridized monolayers using density functional theory (DFT), GW and 
Bethe-Salpeter equation calculations. The CBN optoelectronic properties result from the overall 
monolayer bandstructure, whose quasiparticle states are controlled by the C domain size and he 
at separate energy for C and BN without significant mixing at the band edge, as confirmed by 
the presence of strongly bound bright exciton states localized within the C domains. The resulting 
absorption spectra show two marked peaks whose energy and relative intensity vary with composition 
in agreement with the experiment, with large compensating quasiparticle and excitonic corrections 
compared to DFT calculations. The band gap and the optical absorption are not regulated by the 
monolayer composition as customary for bulk semiconductor alloys and cannot be understood as a 
superposition of the properties of bulk-like C and BN domains as recent experiments suggested. 

PACS numbers: 73.22.-f, 78.67.-n, 71.35.-y 
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Recently, two-dimensional (2D) materials with a finite 
band gap and high carrier mobility have been synthe- 
sized and characterized, chiefly in view of their use in 
optoelectronic devices [il, '^S]. Among these, monolayers 
of hybridized carbon and boron nitride (CBN) present 
highly appealing optoelectronic properties due to their 
tunable optical gap, physically deriving from the large 
band gap difference between pure hexagonal BN (/i-BN, 
energy gap Eg > 5.0 eV) and graphene, a semimetal with 
zero energy gap 0]. The C and BN phases are immisci- 
ble in 2D, leading to phase separation in the sheet with 
formation of distinct C and BN domains jl| . 

The optical absorption spectrum of CBN measured in 
the recent work from Li et al. 1] shows two main absorp- 
tion edges located (using Tauc's extrapolation procedure 
[3|) at around 1.6 eV and 4.5 eV for a sample with 65% 
C content, while for a higher C concentration such ab- 
sorption peaks were shifted to lower energies. These two 
features were attributed to the existence of C and /i-BN 
domains large enough to resume their individual band 
gap identities, rather than behave as a CBN alloy. 

The electronic band gap of the CBN monolayer has 
been calculated for different C and BN domain ge- 
ometries including quantum dot and superlattices that 
approximate the experimental structure, using tight- 
binding and DFT calculations [Hl-flU . 

For quantum dot geometries, the band gap tunabil- 
ity arises due to confinement in both dimensions [sl, but 
DFT level calculations were only applied to limited do- 
main sizes up to wl nm, well below the experimental 



value. CBN 2D superlattices can provide band gap values 
closer to the experimental regime due to lack of confine- 
ment in both dimensions. However, the results appeared 
in the literature up to now have addressed the band gap 
tunability as a function of composition Q rather than 
domain size, or the differences between hydrogenated C 
nanoribbons and C nanoribbons formed within a CBN 
superlattice 0, [12] , without attempting to interpret the 
experimental optical absorption spectrum in Ref. [l| and 
its connection to the calculated bandstructure. 

Here we shed light on the optoelectronic properties of 
CBN monolayers and show that the double optical peak 
measured in [l| is inherently related to the bandstruc- 
ture of the CBN system, regardless of the composition 
and domain size. Such excitation peaks are not due to 
large C and BN domains, but rather to the near-absence 
of mixing of C and BN states at the valence band max- 
imum (VBM) and conduction band minimum (CBM), 
generating two distinct sets of optical transitions: low 
energy transitions between C states near the band gap 
resulting in bound bright excitons localized within the C 
domains, and higher energy transitions from states that 
are mostly BN-like in character and lie deeper withiir the 
conduction and valence bands. 

In addition, we demonstrate that composition is not 
the main variable regulating the CBN electronic struc- 
ture, band gap and optical absorption, in contrast with 
the well-known behavior of bulk semiconductor alloys 
whose energy gap varies continuo usly with the concentra- 
tion of the composing elements [13[. We also find large 
corrections to the DFT band gap using GW calculations 
and strongly bound excitons, both effects contributing to 
an optical absorption spectrum similar to the random- 
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phase approximation (RPA) DFT spectrum due to large 
error compensations in the latter. 

We carried out ab-initio DFT calculations using the 
Quantum Espresso code [3l separately on CBN superlat- 
tices with armchair and zig-zag edges, with a total num- 
ber of C and BN dimer rows of 8 and 16 in distinct sets of 
calculations, so that electronic structure data are avail- 
able for a given composition for two different C and BN 
domain sizes. The Perdew-Burke-Ernzerhof exchange- 
correlation functional is adopted and ultrasoft psue- 
dopotentials [l^ are used to describe the core electrons. 
An orthorhombic unit cell with 32 atoms was used for the 
armchair structures and a hexagonal unit cell with 128 
atoms was used for the zig-zag structures (Fig. [1]), with 
interlayer spacing of 15 A in both cases. A kinetic energy 
cutoff of 35 Ry was used for the plane- wave basis set and 
of 200 Ry for the charge density. Converged Monkhorst- 
Pack fc-grids [17| of 3 x 3 x 1 and 15 x 5 x 1 were used 
for the zig-zag and armchair cases, respectively. 

The GW and Bethe-Salpeter equation (BSE) calcula- 
tions [1^] were performed on three armchair cases [l^l 
with different C concentrations using the Yambo code 
[20| [2l| . Briefly, a plasmon-pole model was adopted 
for the self-energy and cut-off energies of 35 Ry and 5 
Ry were used, respectively, for the exchange and correla- 
tion part of the self-energy; the Coulomb interaction was 
truncated in the direction perpendicular to the sheet to 
avoid spurious interaction with the image system. The 
GW calculations were performed without self-consistency 
in the Green's function and the screened Coulomb inter- 
action (Go Wo approximation scheme). A total number 
of 300 bands (> 150 empty bands) was used, together 
with a converged fc-point grid of 16 x 6 x 1. Both DFT 
and GW levels of theory were employed in combination 
with RPA or BSE calculations to compute the optical 
absorption spectra. 

Fig. [Ha) shows the evolution of the DFT band gap as 
a function of C concentration for armchair superlattices 
with a simulation cell constituted by a total of 8 and 
16 atom rows. For a given composition, the two cases 
with different number of atom rows give different band 
gap values; a same band gap value for the two systems is 
found when comparing superlattices with a same number 
of C rows (and thus a same C domain size). For example, 
the arrows in Fig. [Ija) point at structures with three C 
rows, corresponding to different concentrations for sys- 
tems with a total of 8 and 16 rows in the unit cell, and 
yet showing a same energy gap to within 0.1 eV. This 
behavior is observed for all such pairs of structures with 
a same number of C rows. 

A similar trend is found for the zig-zag case (Fig. [ija)) 
where systems with a given number of C rows, but cor- 
responding to different C concentrations, show a same 
energy gap to within 0.1 eV. As observed in previous 
work the addition of a single C row causes a drop 
in the band gap from 4.6 eV for pure /i-BN to < 2.0 eV 
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FIG. 1. (color online) (a) Kohn-Sham band gap calculated 
using DFT for CBN superlattices with armchair (left) and 
zig-zag (right) edge. Starting from pure BN (zero %C), each 
consecutive point for increasing C concentrations corresponds 
to the addition of a C row to the structure, (b) Unit cells 
employed in the 8-rows calculations for the armchair (left) 
and zig-zag (right) cases. The parts referred to as "one row" 
in the text are shown in dashed boxes. 



within DFT calculations. 

A more detailed analysis of the DFT electronic struc- 
ture is presented in Fig. [21 where the CBN systems are 
labelled Cx(BN)(8_x) for superlattices with x carbon rows 
out of a total of 8 rows. A direct band gap was found for 
all the armchair and zig-zag systems studied, thus justify- 
ing the high optical absorption observed experimentally. 
In the armchair cases shown here, the band gap closes 
progressively for increasing C concentrations but the for- 
mation of the Dirac cone only occurs when 1 — 2 residual 
BN rows are present, as seen by the increasing dispersion 
in the Y — T direction (normal to the domain edge) re- 
sulting in the closure of the gap at the equivalent of the 
K point of the hexagonal lattice . We interpret this 
behavior as a sign of the incipient delocalization of the 
VBM and CBM states when the system is close to being 
a sheet of pure graphene. 

The projected density of states (PDOS) shows that 
for all the compositions studied the states with energy 
near the gap are mainly due to C, while states farther 
in energy from the gap result from a clear hybridization 
of C and BN states. The zig-zag case shows analogous 
behavior (see Supplemental Material ), with subtle 
differences in the band gap closing near the pure C com- 
position that occurs through a band-crossing mechanism 
for the formation of the Dirac cone at the K point, caus- 
ing the slight bump seen in the band gap values in Fig. 
[T^a) at large C fractions. 

The dependence of the band gap on the C domain size, 
rather than on the overall CBN layer composition [20], 
can be regarded as an electronic size effect, whereby the 
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FIG. 2. (color online) PDOS and corresponding bandstruc- 
ture plots for three armchair superlattices. When the number 
of C rows is increased, the band gap decreases continuously. 
However, the Dirac cone typical of graphene only starts form- 
ing for 1 — 2 residual BN rows. The PDOS shows that the 
states close in energy to the gap are mainly due to mixing of 
C valence states with almost no contribution from BN. 



C states near the band gap are confined by an effective 
quasiparticle barrier within the conduction and valence 
bands formed by the BN states, as seen by the analysis 
of the PDOS. Nevertheless, the detail of the edge type 
and binding at the CBN interface is also relevant to de- 
termine the specific value of the gap (and consequently 
of the optical properties), as shown by the different band 
gap values for the zig-zag and armchair cases. 

Next, we analyze the results from the GW and BSE 
calculations. For the three armchair cases studied with 
beyond-DFT methods, we found quasiparticle (GW) 
band gap values significantly higher than the DFT gap 
(Table HI, with corrections in the range of 0.8 — 1.8 eV. 
On the other hand, the optical gap (GW-BSE) values 
were found to be within 0.4 eV of the DFT gaps due to 
large exciton binding energies in the range of 0.7 — 1.5 eV 



(Table HI , inferred from the difference between the quasi- 
particle (GW) and optical (GW-BSE) gaps. 
The DFT-RPA absorption spectra for the armchair (Fig. 



TABLE I. Values of the band gap (in eV) within different 
approximations for three armchair structures studied with 
beyond-DFT methods. The exciton binding energy (in eV) is 
also shown. 



Structure 


DFT 


GW 


GW-hBSE 


Exciton b. e. 


Ci(BN)7 (12.5% C) 


1.39 


3.25 


1.80 


1.45 


C4(BN)4 (50% C) 


0.60 


1.72 


0.59 


1.04 


C7(BN)i (87.5% C) 


0.30 


1.04 


0.32 


0.72 



[3l[a)) and zig-zag (see Supplemental Material) cases show 
two absorption peaks resulting from the DFT bandstruc- 
tures in Fig. [21 with transitions due to the direct band 
gap along the y — F direction (low energy absorption 
onset) and to the direct gap from higher energy states 
formed at the X point. These two features are retained 
in the GW-BSE spectrum (Fig. El^b)) and can ultimately 
be indicated as the reason for the two absorption peaks 
observed experimentally in 

Nevertheless, the GW-BSE spectrum contains infor- 
mation about two important physical effects: the open- 
ing of the quasiparticle gap due to electron-electron in- 
teraction at the GW level of theory, causing by itself an 
almost rigid blue-shift of the absorption spectrum (GW- 
RPA curve in Fig. [U^b)), and the formation of bound 
excitonic states resulting in a red-shift of the spectrum 
back to energies similar to the DFT-RPA approximation 
level, as seen in the GW-BSE spectra in Fig. El^b). This 
explains the excellent qualitative agreement of the DFT- 
RPA optical absorption spectrum seen in Fig. [3][b), due 
to a compensation of large errors within DFT. 

Consistent with the experiments in [ij, for increasing 
C concentrations both absorption peaks are red-shifted 
and their relative intensity varies, thus showing the tun- 
ability of the optical properties of the CBN sheet. By 
comparing with Ref. [l[ where systems with over 65% 
C concentration show an optical gap of 1.5 eV, we in- 
fer that such a low energy absorption edge arises from 
small C domains with at least one dimension of 1 — 2 
nm length. However, we note that the Tauc's extrapola- 
tion procedure Q used in Ref. [l[ can lead to significant 
errors in the estimation of the absorption onset, and sug- 
gest that the experimental spectrum should be regarded 
as formed by the superposition of the spectra shown in 
Fig. [3] weighed for the different C domain sizes present 
in the system, a fact justified by the absence of a well- 
defined low energy peak in the experimental absorption 
data til. 

The exciton wavefunction obtained within the BSE 
framework (Fig. Etc)) is found to be confined by an ef- 
fective quasiparticle gap at the CBN interfaces: when 
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Absorption spectra: £2 (a. u.) vs. Energy (eV) 
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FIG. 3. (color online) (a) Evolution of the DFT-RPA optical absorption (expressed as the imaginary part of the dielectric 
tensor, £2) for increasing sizes of the C domain, for the armchair structures with 8 atom rows shown in Fig. [IJa). The two 
absorption peaks red-shift and vary their relative strength for an increasing C domain size as pointed by the red arrows, (b) 
Comparison of absorption spectra for the armchair structures discussed in Table |I] The imaginary part of the dielectric tensor 
(averaged for the three directions) is shown for different levels of approximation including DFT-RPA, GW-RPA and GW-BSE. 
(c) Exciton spatial distribution (shown in green) for the C4(BN)4 armchair superlattice when the hole is fixed in the position 
represented by the magenta spot in the C domain. The atoms coloring follows the legend in Fig. [Ub). 



the hole is fixed in a given position within the C do- 
main, the electron wavefunction localizes within the same 
C domain, compatible with the presence of VBM and 
CBM states with C character, and BN-like states farther 
in energy from the gap. Such delocalized, yet strongly 
bound bright excitons are a unique feature of nanoscale 
systems that has been predicted previously for ID arm- 
chair graphene nanoribbons [23]. 

In summary, we elucidate the mechanisms underlying 
the band gap formation and optical processes including 
absorption and excitons in the CBN monolayer. We find 
that the energy gap is uniquely regulated by the size 
of the C domains regardless of the overall system com- 
position, a novel mechanism in clear contrast with the 
behavior of bulk semiconductor alloys. The absorption 
spectra of CBN monolayers show two main absorption 
peaks whose strength and position depend on the C do- 
main size, with transitions between the VBM and CBM 
states contributing to the low energy absorption onset 
and to the formation of strongly bound excitons within 
the C domains. Large quasiparticle and excitonic correc- 
tions show that DFT calculations are inadequate to pre- 



dict the experimental behavior of 2D CBN alloys. The 
unique dependence of the band gap and absorption spec- 
trum on domain size rather than concentration opens 
new possibilities for band gap engineering and in general 
for optoelectronic and photovoltaic applications distinct 
from those existing for /i-BN and graphene. Similar large 
excitonic effects have been predicted recently for M0S2 
monolayers [l^, suggesting that this new family of 2D 
nanomaterials could constitute a novel playground for 
optoelectronic and excitonic devices. 
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* E-mail: jcg@mit.edu' 

[1] L. Ci, L. Song, C. Jin, D. Jariwala, D. Wu, Y. Li, A. Sri- 
vastava, Z. F. Wang, K. Storr, L. Balicas, F. Liu, and 
P. M. Ajayan, Nature Mater. 9, 430 (2010). 

[2] B. Radisavljevic, A. Radenovic, J. Brivio, V. Giacometti, 
and A. Kis, Nature Nanotech. 6, 147 (2011). 



5 



[3] 
[4] 

[5] 
[6] 

[7] 

[8] 

[9 

[lo: 
[11 

[12 

[13 
[14 



A. Geim, Science 323, 1530 (2009). 

J. Tauc, R. Grigorovici, and A. Vancu, Phys. Stat. Sol. 
15, 627 (1966). 

G. Seol and J. Guo, Appl. Phys. Lett. 98, 143107 (2011). 
K.-T. Lam, Y. Lu, Y. P. Feng, and G. Liang, Appl. Phys. 
Lett. 98, 022101 (2011). 

A. K. Manna and S. K. Pati, J. Phys. Chem. C 115, 
10842 (2011). 

J. Li and V. B. Shenoy, Appl. Phys. Lett. 98, 013105 
(2011). 

A. Y. Liu, R. M. Wentzcovitch, and M. L. Cohen, Phys. 
Rev. B 39, 1760 (1989). 

J. R. Martins and H. Chacham, ACS Nano 5, 385 (2011). 
X. Blase, J.-C. Charher, A. D. Vita, and R. Car, Appl. 
Phys. A 68, 293 (1999). 

Y.-W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. 
Lett. 97, 216803 (2006). 
M. Jaros, Rep. Prog. Phys. 48, 1091 (1985). 
P. Giannozzi, S. Baroni, N. Bonini, M. 
R. Car, C. Cavazzoni, D. Ceresoh, G. L. 
M. Cococcioni, I. Dabo, A. Dal Corso, S. 
coh, S. Fabris, G. Fratesi, R. Gebauer, U. Ger- 
stmann, C. Gougoussis, A. Kokalj, M. Lazzeri, 
L. Martin-Samos, N. Marzari, F. Mauri, R. Maz- 
zarello, S. Paohni, A. Pasquarello, L. Paulatto, 
C. Sbraccia, S. Scandolo, G. Sclauzero, A. P. Seitso- 
nen, A. Smogunov, P. Umari, and R. M. Wentzcovitch, 
'Journal of Physics: Condensed Matter 21, 395502 (19pp) 



Calandra, 
Chiarotti, 
de Giron- 



[15] J. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 

77, 3865 (1996). 
[16] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 
[17] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 

(1976). 

[18] G. Onida, L. Reining, and A. Rubio, Rev. Mod. Phys. 
74, 601 (2002). 

[19] We calculated the formation energy of the armchair and 
zig-zag type edges and inferred a ratio of 2.5 between 



the thermal equilibrium length of the armchair and zig- 
zag domain edges at intermediate compositions and room 
temperature. This justifies our more accurate study of the 
optical properties and band gap for the armchair edge 
type as it is expected to contribute more significantly to 
the experimental results observed in 

[20] A. Marini, C. Hogan, M. Gruning, and D. Varsano, 
Comp. Phys. Comm. 180, 1392 (2009). 

[21] For the calculations performed with the Yambo code, 
the ground state Kohn-Sham wavefunctions and eigen- 
values were obtained using an LDA exchange-correlation 
functional (Phys. Rev. B 23, 5048, 1981) with TrouUier- 
Martins norm-conserving pseudopotentials (Phys. Rev. B 
43, 1993, 1991) as implemented in the Quantum Espresso 
code. A kinetic energy cutoff of 55 Ry was used for the 
wavefunction. These calculations yielded identical results 
to the ones using the PBE exchange-correlation func- 
tional with ultrasoft pseudopotentials, and were carried 
out with the only scope of obtaining an input compatible 
with the Yambo code. 

[22] See Ref. Q for the nomenclature of the Brillouin zone 
for the orthorhombic cell used here. The K point of the 
hexagonal lattice, where the graphene band gap is ex- 
pected to close at high C concentration, falls between 
the r and Y point within this scheme. 

[23] See Supplemental Material at URL for the bandstructure 
and DFT-RPA optical absorption spectra for the zig-zag 
'(2009)1 CBN system, and for supplemental evidence of the de- 
pendence of the band gap on C domain size rather than 
composition. 

[24] X. Zhu and H. Su, J. Phys. Chem. C 114, 17257 (2010). 
[25] T. Olsen, K. W. Jacobsen, and K. S. Thygesen, e-print 
arXiv: 1107.0600 (2011). 



